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We present an analysis of the magnetic field distribution in the Abrikosov lattice of high-K su- 
perconductors with singlet pairing in the case where the critical field is mainly determined by the 
Pauli limit and the superfluid currents partly come from the paramagnetic interaction of electron 
spins with the local magnetic field. The derivation is performed in the frame of the generalized 
Clem variational method which is valid not too close to the critical field and furthermore with the 
Abrikosov type theory in the vicinity of it. The found vortex lattice form factor increases with in- 
creasing field and then falls down at approach of the upper critical field where the superconducting 
state is suppressed. 
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I. INTRODUCTION 

Recent small angle neutron scattering experiments per- 
formed on the heavy-fermion superconductor CeCoIn 5 
have revealed an unexpected behavior of the vortex lat- 
tice (VL) form factor-^] defined as the Fourier transform 
of the local magnetic field in the vortex lattice. The VL 
form factor of the type II superconductors is usually a de- 
creasing function of the magnetic fieldP On the contrary, 
the VL form factor in CeCoIns was found to increase 
with increasing magnetic field and the n to fall down at 
the approach of upper critical fieldP^ 

CeCoIns is a tetragonal, d-wave pairing superconduc- 
tor with a large Ginzburg-Landau (GL) parameter (i.e. 
k = A/£ 3> 1, where A and £ are the two characteristic 
lengths of the GL theory) , and with the highest critical 
temperature (T c = 2.3 K) among all the heavy fermion 
compounds.™^ It has already generated great interest 
caused by the signs of the existence of the FFLO phase 
for a magnetic field parallel to the ab plane (and possi- 
bly to the c-axis)P and by the presence of a temperature 
interval < T < Tq where the superconducting/normal 
phase transition is of the first order The tempera- 
ture To was found equal to 0.3T C for the field orientation 
parallel to the tetragonal axisP 

The distance between the vortices in the mixed state of 
type II superconductors decreases with increasing mag- 
netic field. It diminishes with respect to field penetra- 
tion depth. As a result the field distribution is smoothed 
towards H c i- Hence the form factor characterizing the 
sharpness of local field modulation is also inevitably di- 
minished with increasing field. The unusual field depen- 
dence of the form factor in CeCoIns points out that some 
other mechanism leads to the opposite process, namely to 
the increasing of fast field modulation with field growth 
(Fig. [T]). The anomalous magnetic field dependence of 
the flux lattice form factor observed by neutron scatter- 
ing has been attribut ed to t he large paramagnetic depair- 
ing effect in CeCoIn 5 P±2] 

The relative weight of the orbital and the Zeeman 
mechanisms for suppression of the superconducting state 



by a magnetic field is quantified by the Maki parame- 
ter a m — V2H C 2o/H p , where H&o = <Po/2tt^q is the 
orbital critical field while H p = Ao/v^A* is the Pauli 
limiting field. (j> ~ 2.07 x 1CT 7 G ■ cm? is the flux quan- 
tum. The Maki parameter is expressed through the Fermi 
velocity vp = kp/m* and the critical temperature as 
OiM ~ T c /mv F , hence it takes larger values in heavy 
fermionic compounds where effective mass of charge car- 
riers is much larger than the bare electron mass m* ^> m. 
Unlike the majority of superconductors, the Pauli limit- 
ing field in CeCoIns is smaller 9 than the orbital critical 
field by a factor of > 3. Hence, the Zeeman interaction 
plays an important role in the mixed state field and cur- 
rent distributions. 

The role of the paramagnetic mechanism has been 
investigated using numerical proce ssing of the quasi- 
classical Eilenberger equations! 11 * 12 ! In the first papeJ^ 
the authors could not reproduce the form factor increas- 
ing behavior towards H C 2 observed in CeCoIns. Then 
in the subsequent calculation^ performed with an even 
stronger paramagnetic contribution they have obtained 
the desirable correspondence. The developed numerical 
procedure has been performed at temperature T = 0.2T C 
where the phase transition from the normal to the super- 
conducting state is of the first order. In this region of 
the phase diagram an analytical treatment of the prob- 
lem is impossible. However, at temperatures higher than 
To and at large Maki parameters one can develop an an- 
alytic theory giving the possibility to compare directly 
the roles of the orbital and the Zeeman contributions in 
the form factor field dependence. 

Our approach is based on the ClerrP elegant method 
which approximates the GL order parameter by a trial 
function and allows an analytic solution of the second GL 
equation for the local magnetic field. Here we present an 
analytic derivation of the magnetic field distribution and 
VL form factor taking the paramagnetic effects into ac- 
count. For this pu rpose we shall use the GL theory devel- 
oped in the paper^IHl where the diamagnetic superfluid 
currents are determined not only by the orbital effects 
but also by the Zeeman interaction of the electron spins 
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with the local magnetic field. The Clem-type analysis 
which makes use of the isolated vortices approximation 
is valid not too close to the critical field. Near the critical 
field it should be completed by the magnetic field distri- 
bution based on the Abrikosov type distribution of the 
order parameter. 

By direct solution of the Maxwell equation including 
the currents due to electron spin interaction with mag- 
netic field we demonstrate analytically that the VL form 
factor growth is related to sharp changes of the local mag- 
netic field concentrated in the cores of vortices that is on 
a distance of the order of the coherence length around 
the vortex axis. The found form factor decreases with in- 
creasing magnetic field in the high temperature - low field 
region of the phase diagram, while at lower temperatures 
it turns to the behavior increasing with increasing field 
and then decreasing to zero at the approach of th e critical 
field. In accordance with numerical studies^^ we have 
found that the non-monotonic behavior reveals itself in 
case of strong enough paramagnetic contribution. The 
obtained results qualitatively describe the magnetic field 
dependence of the vortex lattice form factor in any high-K 
superconductor with singlet pairing and large Maki pa- 
rameter. At the same time we do not pretend on a quan- 
titative correspondence of our findings with the CeCoIns 
form factor dependence reported in the paperJ^. In 
absence of knowledge of real band structure, the Fermi 
surface shape, pairing mechanism, probable field depen- 
dence of the effective mass, the value and the angular de- 
pendence of the (/-factor, any t heory including the devel- 
oped numerical procedure^^ cannot pretend on a quan- 
titative correspondence with experiment. Therefore, our 
goal is to develop an analytical model which gives us a 
clear picture of the field distribution in the Abrikosov lat- 
tice of the superconductors with a large Maki parameter 
value. 

The paper is organized as follows. In the Second Sec- 
tion we briefly repeat the Clem results and discuss the 
limitation of their applicability. Then, in the third Sec- 
tion the corresponding theory taking into account the 
Zeeman interaction is presented. The behavior of the 
form factor in the GL region near the critical field is de- 
scribed in Section IV. The Concluding remarks are for- 
mulated in Section V. 



II. ORBITAL FORM FACTOR 

Regarding the form factor field dependence, the par- 
ticular symmetry of the flux line lattice is unimportant. 
For simplicity we shall consider a square vortex lattice 
with axis spacing a = y/fa/B formed in a type II super- 
conductor under magnetic field directed along the z-axis. 
The magnetic induction B — h is determined as the spa- 
tial average of the local magnetic field h = hz = V x A. 
The form factor F mn is determined by the Fourier trans- 



form of the magnetic field 



h(q m „) = / d 2 r h(v)e 



-»q„ 



(1) 



where q m „ = — (mx + ny) are the vectors of the re- 
ciprocal lattice. For the external field not too close to 
H C 2, where the distance between the vortices is larger 
than the core radius, the local magnetic field repre- 
sents the sum of the magnetic fields of separate vortices 
h(r) = J^i h v {v — rj). Thus, the form factor is propor- 
tional to the Fourier transform of magnetic field around 
one vortex 



B 

00 



(2) 



This quantity is related to the intensity of Bragg peaks 
observed in small angle neutron scattering experiments. 
Throughout the article we shall consider the form factor 
Fqi that corresponds to the indices (0,1). We designate 
it as F. 

The GL theory for the VL form factor, valid in the 
limit k 1 was developed by J. ClemP Starting from 
the general form of the order parameter for an isolated 
vortex 



A(r) = A 0o /(r)( 



(3) 



(ip is the azimuthal angle in the plane perpendicular to 
the vortex axis), he proposed to model f(r) by the trial 
function 



r 



(4) 



with R = yV 2 + The variational parameter was 
constrained to minimize the vortex total energy and was 
found to be in the large k limit = \J2 £, where £ is the 
coherence length. Clem has calculated the field distri- 
bution due to the orbital current and obtained the form 
factor 



F nr h — B 



(5) 



where Q = \/q 2 + A~ 2 , q = 2ny/B/(j) , Ki(x) is the 
modified Bessel function of first orderf^ and A is the 
London penetration depth. One can write an approxi- 
mate form of it in the conditions k> 1 and <; > A -1 , 



orb 



7 q£ v Ki(q£ v ). 



(6) 



The found form factor slowly decreases with magnetic 
field. While the product q£ v is smaller than 1, that is 
not too close to the upper critical field, the form factor 
linearlv decreases with field. The formal application of 
Eq. (|6j) up to H fa H c2 , where q£ v ps \/2~Tt, gives the 
exponential decrease of the form factor. In fact in the 
vicinity of H C 2 the approximation of independent vortices 
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does not work. The proper calculation should be done 
using the Abrikoso\051 form of the order parameter (see 
Section IV) that leads to the vanishing of the form factor 
linearly with (H C 2 — H). 

Equation ([5| is valid in the region of validity of the GL 
theory. The latter does work near the critical tempera- 
ture (T c — T)/T c <C 1 but not at low temperatures where 
the gradient expansion of the free energy is inapplicable 
and one needs to take into account the higher order gra- 
dient terms. Despite of this fact equation ^ is widely 
used in the discussions of the form factor field depen- 
dence at low temperatures T<T C (see for instance Ref. 
3, 17). The higher order gradient terms change the con- 
crete coordinate dependence of the order parameter near 
the vortex axis but details of this dependence have no 
significant influence on the field distribution in the vor- 
tex lattice. Hence, in the strong type II superconductors 
even at low temperatures Eq. ([5| still gives the quali- 
tatively correct description of the form factor magnetic 
field dependence for the fields not too close to H c i. 

III. THE ZEEMAN CONTRIBUTION 

The form factor given by equations (J5j , is reliable if 
one neglects the paramagnetic interaction of the electron 
spins with the magnetic field. The latter leads to two ex- 
tra features that are existent in the case of a large enough 
Maki parameter. First, the two characteristic lengths £ 
and A in the above expression prove to be slightly mag- 
netic field dependent. Second, a new mechanism originat- 
ing from the Zeeman interaction gives rise to a contribu- 
tion to the diamagnetic screening-^ in the high magnetic 
field region of the phase diagram. To find it we start with 
the Ginzburg-Landau formulation including the param- 
agnetic effects. 

The superconductor CeCoIn 5 has pairing symmetry 
d x 2_ y 2, ls with order parameter 

A k (r) = V(k)A(r), ^(k) = v^cos^). (7) 

For s-wave superconducting state ip(k) = 1. The free 
energy of the system is given by the GL functional 

T = J d 2 r ( ^ + a|A| 2 + e{h z - B)|A| 2 

+/ 3|A| 4 + 7 |DA| 2 ), (8) 

where D = — iV + 2eA is the gauge- invariant gradient 
(from here we put H = c = 1), h = rot A is the local 
internal magnetic field, and the coefficients in the func- 
tional depend on both temperature T and induction B 
determined by the spati al average h = B = Bz. In the 
clean limit they ar d 13 * 14 l 

a = N (]n(T/T c )+Ut^{w)-^(l/2)), 

da N Q fi . 
e = —— = — — isrnw (in), 
dB 2ttT v h 



6 = -^^(l^(k)| 4 )5Re*( 2 )H, 

7 = -g^a<|^(t)| a ^>»e* (a) H, 

where ty(w) is the digamma function, ^^(w) are its 
derivatives called by the polygamma functions,^ and 

1 iuB 
w = . 

2 2ttT 

We shall consider s- and d-wave superconducting states 
in a quasi-two-dimensional crystal with a nearly cylin- 
drical Fermi surface. For the <i-wave order parame- 
ter given above the averages over the Fermi surface are 
(|^(k)| 4 ) = 3/2 and (|V>(k)| 2 fc 2 )=l/2 while for s-wave 
superconductivity and the same Fermi surface the corre- 
sponding averages are 1 and 1/2. 

In the paramagnetic limit, when the orbital effect is 
neglected the transition from the normal to a supercon- 
ducting state takes place at the critical field B C (T) de- 
fined by equation a(T, B) = 0. Along this transition line, 
the coefficients /3(T,B) and "f(T, B), which are positive 
near T c , become negative at T < T* ~ 0.56T C . This de- 
fines the tricritical point (T*,B*) of the phase diagram 
with B* = B C (T*) ~ 1.07T c //i. At the tricritical point, 
the sign change of the coefficient 7 signals an instabil- 
ity toward the FFLO state with spatial modulation of 
the order parameter A, while the sign change of the co- 
efficient j3 signals a change of the order of the normal 
to superconducting phase transition. A more elaborate 
treatment including the orbital effects and the higher or- 
der terms in the GL functionate results in the following 
effects: (i) the upper critical field is slightly reduced by 
value of the order B C (T) / 'au\ (h) the temperature where 
the change of the order of the transition occurs and the 
one where the FFLO state arises are decreased by values 
of the order T c /aM with respect to T* . Below we con- 
sider only the temperatures above T* where the j3 and 
7 coefficients are positive and it is not necessary to take 
the higher order gradient terms into account. 

In the case of a large Maki parameter the GL expansion 
of the free energy in powers of the order parameter and 
its gradients is valid near the critical field which i s main ly 
determined by the paramagnetic depairing effect At 
smaller fields strictly speaking the higher order gradient 
terms should be included. The situation is similar to the 
application of the Clem formula at low temperatures dis- 
cussed in the previous Section. The higher order gradient 
terms change the concrete coordinate dependence of the 
order parameter near the vortex axis but details of this 
dependence have no significant influence on the field dis- 
tribution in the vortex lattice. Hence, in the strong type 
II superconductors even for the fields noticeably smaller 
than the critical field the calculation with a GL functional 
containing just the second order gradient term gives the 
qualitatively correct description of the form factor mag- 
netic field dependence.^ 
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Following Clem's procedure we consider an isolated 
vortex with order parameter given by Eqs. @ and 
amplitude Aoo = y/\a\/2/3. The field distribution around 
a single vortex is determined by the Maxwell equation de- 
rived from the stationary condition of the GL functional 
with respect to the vector potential 



An 

The density of current 



-V x h„ 



: 3orb + Jz 



(9) 



(10) 



consists of two parts originating from two different terms 
in the GL functional. The orbital density of current is 



jorb = -8e 2 7 [A v (r) 
while the Zeeman current!^ is 
jz 



2nr 



|A| 2 ^, 



dr 



(11) 



(12) 



Here the vector potential has the form A„(r) = A v (r)(p. 
Hence, with the help h„ = V x A„, we come to the 
equation that determines the vector potential A v (r) 



dr r dr v 



t 
A 2 



A, 



df 2 

, v , AneAl^-, (13) 
2nX 2 r dr 



M 2 



where A = y /3/167re 2 7|a| is the penetration depth. 
Let us introduce the auxiliary function 



v s (r) 



9o 
2nr 



Mr) 



playing the role of supernuid velocitj^ and substitute 
this into Eq. (13). We obtain the differential equation 



with an inhomogeneous term of Zeeman origin 

dr { rdr [ s)) A 2 Vs ~ 47r£A °° dr ' 
The general solution for this equation 

v s {r) = vi(r) + v h s {r) 



(14) 



(15) 



consists of the sum of particular solution of the inhomo- 
geneous equation ( 14 1 and a solution of corresponding 



homogeneous equation. The former is given by 
R 



where 
C(z) 



vl(r) = - 
A 



-K X {R/X)C{R/X), 



(16) 



dz 



KM 



dz, (17) 



chosen such that t>*(0) = and v\(oo) — 0. The latter 
condition is assured by letting the constant be zero in the 
primitive j z dz of the function Kl i z ^ ■ 



The falling to zero at r — >• oo solution of the homoge- 
neous equation 



v h s (r) = 



RKi(R/X) 



meets the requirement that the vector potential 



Mr) 



2irr 



v h s {r) - v\{r) 



(18) 



(19) 



vanishes on the vortex axis. 

One can divide the total vector potential by its orbital 
part and its Zeeman part 



Mr) =A orb (r)+A z {r), 
where the orbital part 



(20) 



-WO (2D 
2nr 2nr \ £vKl\£v/X)/ 



is the solution of the equation ( 13 ) in the absence of the 



last term of the Zeeman origin that was found in Ref. 4. 
The corresponding magnetic held h or b = h or f,z is 



foorb — 



0o K (R/X) 
2wX^K 1 ^ v /Xy 



(22) 



and the form factor is determined by Eq. |5]). 

The Zeeman part of the vector potential is given by 

A z (r) = -vl(r) = -Ki(R/\)C{R/\). (23) 
r 

The corresponding magnetic held hz = hzz reads 

h z {r) = - [-K (R/X)C{R/X) + K 1 {R/X)C'{R/X)] . 

The numerically found magnetic held coordinate depen- 
dence is shown in Fig. 1. 

It is worth noting that the Zeeman term does not spoil 
the basic properties of the Abrikosov vortex. Namely the 
total magnetic flux through the surface perpendicular to 
the vortex axis is equal to the flux quantum 



2tt / (h or b(r) + h z {r))rdr = 
Jo 

2ir lim r{A orb {r) + A z {r)) = O (25) 

r— loo 

To prove this property, one must take into account the 
asymptotic behavior at large distances of functions 



Ki{z) 



z > 1, 



and 



C(z) 



A 



2-5/2' 



z > 1. 



(26) 



(27) 



5 



The Zeeman part of the vector potential Az(r) is deter- 
mined by the product of these functions and diminishes 
at r 3> A as oc 1/r 3 . 

On the other hand at small distances these functions 
behave as 



Ki{z) 



z < 1, 



and 



C7(z)«^^ln^, 

A fu 

For the superfluid velocity given by 



z < 1. 



(28) 



(29) 



«s(r) 



27rr VC^ifeA) O 
at r -C A we obtain 



1 



In 



7? 



(30) 
(31) 
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The magnitude of dimensionless combination 

is estimated as follows. Using the relations Aj^^ 
and e = N fj^m^' (w)/2nT rj 7C(3)m*fc F ^ 2 B/(27r 2 ) 2 T 2 
(we remind that u> = 



we come to 



^ 2 e&Ul 7C(3) . 
r = — ^-fc^re — — <1. (32 




FIG. 1: Plot of the radial distribution of the magnetic field 
around a single vortex at induction B = LOT and temper- 
ature T = 1.05T* calculated using parameters pointed out 
in text. The blue curve represents the orbital part of mag- 
netic field (see the text) while the green one corresponds to 
the total field. The field hz(r) has significant variation near 
the vortex core. While at large distances of the order of inter 
vortex spacing (not shown here) it starts to be negative, such 
that the total flux around single vortex is still equal to flux 
quantum (see text). 



Here Uf is the Fermi momentum, r e = e 2 /mc 2 is the clas- 
sical radius of the electron, the product of these values 
is kpr e f=a 10 -5 . The ratio of the Fermi energy in heavy 
fermionic compounds to the temperature Ef/T taken at 
temperatures T ~ T c /2 is not larger than 10 2 -=- 10 3 . The 
ratio [J.B/T rs \xB /T c is smaller than 1 even in the com- 
pounds with heavy electron mass. Thus, inequality (32) 
is certainly carried out. 

Hence, we come to the usual expression for the super- 
fluid velocity 



v s (r) 



9o 
2nr 



(33) 



which is valid at r -C A. Furthermore, as far as we are 
working at 5 > H cl the Zeeman term e(h v z — B)A 
yields just negligibly small correction to the regular term 
a(B)A in the GL equation. It means that the solution 
of the GL equation for the isolated vortex has the coor- 
dinate dependence 



|A(r)| 
A™ 



1 



r > S, 



as it is in the absence of the Zeeman interaction.^ Here, 
£ = \Jl/\ a \ i s the coherence length. This justifies the 
variational approach making use the order parameter 
given by Eqs. (|3j), Q. 

The total field amplitude h or b{r) + hz(r) differs sig- 
nificantly from its orbital part h or b(r) only at the dis- 
tances of the order of £ v from the vortex axis (see Fig. 
1). Hence, the energy of a single vortex is dominated by 
the orbital field h or b{r). Therefore the minimization of 
the energy of a single vortex gives the same variational 
parameter £ v = \/2£ as exposed in the Clem paper P The 
lower critical field keeps the usual value determined with 

logarithmic accuracy as H c \ 4 ^° 2 In re. 

In the high-re limit, we can find a simpler expression 
to the Zeeman internal field by looking at its behavior at 
distances r « i By using Eqs. (28) and (29) we find 
the dominating term in Eq. ( 24 ) 



4-TreA;: 



R 2 ' 



(34) 



We see that hz(r) has significant variations near the 
core of a vortex, the characteristic length associated to 
it being of the order of £ v (see Fig.l). Also the local 
field amplitude increases with increasing field because 
e is proportional to B. The formula (34) may be de- 
rived by considering the equation Vxh^ = 47rj2, that 
is valid in the absence of the orbital current. From 
Eq. (24) one can find the correction to Eq. (34), 
Sh z = -iireA^KoiR/X) hx(R/£ v )/ K 2 . 

The fast variations of local magnetic field given by the 
Zeeman interaction in the vicinity of the vortex axis dras- 
tically changes the form factor field dependence. One can 
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calculate the Fourier transform of the the Zeeman inter- 
nal field around a single vortex 



h z (q) = / & z vh z (r)e 



= 2tt / rdrJo(qr)hz(r) 
Jo 

(35) 

using Eq. (34 1 valid for r <C A. Indeed, one can prove 



that the correction to the Fourier transform brought by 
the region r > A is negligibly small if g _1 -C A (that is 
when the distance between the vortices is much smaller 
than the penetration depth). Furthermore the accuracy 
of the approximation ( 34 ) was validated by numerical 
integration of the full expression ( 24 ) . Therefore in what 



follows we can work with expression ( 34 ) instead of the 
more cumbersome Eq. ( 24 1 , the developed approach is 



applicable to the description of the VL under the fields 
which are much larger than the lower critical field and 
not too close to the upper critical field. 

By using the coefficients of the Ginzburg-Landau de- 
composition we obtain the full expression for the 
coherence length 



v 2 F (\^(k)\ 2 kl)^e^ 2 \w) 



8(2ttT) 2 [ ln(T/T c ) + ^(w) - (1/2)] 



(36) 



For the fields far enough from the critical field this for- 
mula gives a slight decrease of the coherence length with 
increasing field before yielding a strong coherence length 
increase below the critical line. 

The contribution to the form factor that originates 
from the interaction of the electron spins with the lo- 
cal magnetic field for an array of B/<fio vortices per cm 2 
is 



8tt 2 £A 2 £ 2 

Fz = -h z (q) = £ A °°^ BK (qt v ). (37) 



D 



After substitution of the explicit expressions of all the 
values it is 

4irN vl uB (W\i)\ 2 k 2 ) ., s 

Fz = ,° ET )' , , 5m$' (w )K >/2gfl, 38 

<t>oT (|^(k)| 4 ) 



where 



we 



used the 



A 2 e 



^(|^(k)| 2 fc 2 )/2(|^(k)| 4 ). 



relation 

FMV-'WI "> x ii*\\v\*>)\ /• Tne g a P averages corre- 
sponding to the s- and d-wave cases around a cylindrical 
Fermi surface are pointed out in the beginning of the 
Section and show that the form factor is bigger in the 
s-wave case than in the d-wave case by a factor of 3/2. 

Let us also rewrite the form factor orbital contribution 
given by Eq. ^ with the coefficients of the model 



Forb 



4irv 2 F (Wk)| 2 fc 2 

00 (iv(k)| 4 ) 



|a|v^g^i(>/2flf0- (39) 



For [iB < 2nT we have 3m$'(w) ~ 7 C{Z) ^B / itT and 
\a\ ~ N (HT c /T)-7((3)(liB/2ttT) 2 ). The ratio of the 
form factors d3Sb and ((39l becomes 



F z 
F, 



K{Z)V 2 B 2 



K (V2q0 



ttT 2 V^qCK^qOiHTc/T) - 7((3)(pB /2ttT) 2 ] ' 

here C( x ) is the Riemann zeta function. For \/2q£, ~ 1 we 
observe that the Zeeman part of the form factor prevails 
over its orbital part in the phase diagram region where 
fiB ~ T. It becomes clear why such a type of behav- 
ior occurs only in the superconductors with dominant 
paramagnetic depairing effect. Indeed at temperatures 
T ~ T c /2 one can estimate nB/T ~ B/H p . For a su- 
perconductor characterized by a small Maki parameter 
the magnetic field B does not exceed the orbital critical 
field B ~ H c2 q, therefore the effect of paramagnetism 
becomes unobservable. The ratio of the form factors de- 
pends from Fermi velocity only through the coherence 
length. The dominant role played by the Zeeman part of 
the form reveals itself at small enough coherence length 
which corresponds to the sufficiently small Fermi veloc- 
ity. 

We have plotted in Fig. [2] the total form factor 



F = F, 



orb 



F z 



(41) 



where F or b and F z are given by Eqs. ((HJ and @ cor- 
respondingly. In numerical calculations we assumed the 
values fi = gfig/2 = \lb for the electron magnetic mo- 
ment in the material, and vp = 5 x 10 5 cm/ s for the Fermi 



First order line 



Critical end point (T+.H*) 




orb 



(40) 



FIG. 2: (Above) CeCoIns phase diagram for B||c-axis. The 
color lines represent the temperatures at which we applied the 
model. (Below) Variations of the squared form factor F 2 at 
different temperatures including both the orbital and Zeeman 
contributions. The dashed lines represent the variations of the 
orbital part only. 
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velocity inside the superconducting phase. A value for vf 
slightly bigger was given in Ref. 10 as a result of mea- 
surements of the upper critical field H c2 near T c . The 
2-D density of states on the Fermi surface is given by 
iVo = m* /2tt£ c , where we considered to* — 100 m e for 
the electron effective mass, and £ c = 7.6 x 10 -8 cm is 
the lattice c-axis spacing. In the form factor variations, 
there is first a domination of the orbital part in the low 
magnetic field region (Fz vanishes at B = 0). We observe 
next a crossover to a region where the paramagnetic term 
is dominant. The regime where ^ goes exactly against 
( 38 ) is likely to explain the observed constant logarithm 
of the squared form factor^ in the interval B = 0.5 — 2T. 
In addition, these anomalous form factor variations were 
observed in the experiments realized on the s-wave su- 
perconductor TmN^I^CP^ 



IV. VICINITY OF PHASE TRANSITION 

The observed form factor^^l falls towards zero near the 
phase transition to the normal state. Its field dependence 
in this region is out of applicability of the derivation that 
makes use of the assumption of isolated noninteracting 
vortices. The form factor decrease at the approach of 
the upper critical field can be found in the frame of the 
GL theory valid in the vicinity of the transition at tem- 
peratures above the tricritical point. 

In vicinity of the phase transition the local magnetic 
fieloP^ deviates from its average value B = h(r) by 



/i 1 (a:,y) = -47r e (|A( a :,y)|' J -|A|a). 



(42) 



The order parameter solves the linearized Ginzburg- 
Landau equation. For a square vortex lattice with period 
a = y/jr/eB it is 



A(x,y) = C (-l) n e 27rmy/a i 



^-(x-na+a/2) 2 



n— — oo 



(43) 

so that \A(ma, na)\ = (m and n are integer). Here, 
\b = l/y2eB is the magnetic length. Multiplying the 
expression by its complex conjugate we have 



+oo 



\A(x,y)\ 2 = C 2 (-l)" +m e" 



^-(x-na+a/2) 2 



771 .71 — — OO 



— — Kj- (x— mo+a/2) 2 o ■/ \ / 
xe 2A| V t 1 } e 2ni{n-m)y/a 



(44) 



Now, putting the dummy index n' = n — m, and using 
Poisson's summation formula^ in the form 



+oo 1 +oo „ 

E = , E fP?-V 



■niraxj a 



where 



a * — ' a 

m— — oo 



f(q)= dxf(x)e-^ x 



is the Fourier transform of f(x), we get the quantity of 
interest 



|A(x,y)| 5 



C 2 



E (-D 



m-\-n-\-mn 



m.n— — oo 



v/2. 

Its average value is 



r (m +n ) 27ri(mx-\-ny) / a 



(45) 



CP 
V2' 



(46) 



Taking into account Eqs. (451 and (461, the Fourier 



transform of the magnetic field ( 42 ) yields the form factor 



corresponding to the Bragg peak with indices (m,n) 



F mn = 47r£|A| 2 (-l) m+n+mll e-^ m +n \ (47) 

The field dependence of |Ap at B ~ H c2 (T) is knownP 
We write it in the limit of large GL parameter k 3> 1 



IAI 2 



e{H c2 - B) 



(48) 



Here e = e{H c2 ) and j3 A = |A| 4 /|A| 2 is the Abrikosov 
parameter. The critical field 

H c2 (T) = B C (T)(1 - 2e 7 / £ ) 

is somewhat lower than B C (T) determined by the equa- 
tion a(B, T) = 0.121 Finally we obtain for the form factor 
near the phase transition line 



F=\F 01 \ 



27re-2£ 2 (ff c2 - B) 



(49) 



When the transition becomes first order but near the 
tricritical point the form factor keeps the form of the 



one calculated and given by Eq. (47). However, the 



main difference here is that the order parameter takes on 
a finite value at the transition. Hence, the form factor 
discontinuously falls down to zero at the transition to the 
normal state. 



V. CONCLUSION 

Making use of the generalized Clem approach we have 
calculated the magnetic field dependence of the vortex 
lattice form factor. The interaction of the electron spins 
with the space inhomogeneous magnetic field existing 
inside the superconductor in the mixed state produces 
diamagnetic currents of pure paramagnetic origin.^ The 
corresponding field concentrated around the vortex cores 
yields a dominant contribution to the expression of the 
vortex lattice form factor at high magnetic field in the 
superconductors with a small enough coherence length 
(that was related here to the Fermi velocity). Finally, 
our analysis showed that the form factor behaves as the 



square of the order parameter in the region close to the 
upper critical field and therefore it falls down to zero at 
the transition to normal state. The fall must be discon- 
tinuous when the transition becomes first order. These 
results give a physical mechanism for the occurrence of 
anomalous VL form factor variations and possibly ac- 
count for the measurements that were recently made on 
CeCoIn 5 . 
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